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ABSTRACT 

The  concept  of  a circular  design  is  defined  and  when  proper 
balance  for  various  effects  Is  assumed.  Its  universal  optimality 
Is  proved  over  the  class  of  all  designs  with  the  same  set  of  para- 
meters. Such  designs  are  shown  to  minimize  the  variance  of  the  best 
linear  unbiased  estimators  of  contrasts  of  residual  and  direct 
effects  over  the  class  of  equlrepllcated  designs.  All  models 
assume  first  order  residual  effects  and  are  of  a circular  nature. 
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1.  intrpdttctipp. 

An  experiment  in  which  a unit  Is  exposed  repeatedly  to  a 
sequence  of  treatments  Is  called  a repeated  measurements  design 
(RMD).  The  experiment  Is  based  on  t treatments « n experimental 
units  and  p periods,  each  unit  being  given  one  treatment  during 
each  period.  A repeated  measurements  design  can  therefore  be 
thought  of  as  a pxn  array,  d,  with  treatments  as  entries.  We 
label  the  treatment  l,2,...,t.  The  entry  In  d Is 

referred  to  by  d(l,J),  l.e.,  d(l,J)  Is  the  treatment  administered 
during  the  1 period  to  the  J experimental  unit.  The  res- 
ponse models  that  we  assume  here  differ  somewhat  from  the  usual 
models  with  residual  effects  (see  Hedayat  and  Afsarinejad  (1973) 
and  (1978),  Magda  and  Hedayat  (1979)*  or  Cheng  and  Wu  (1979)  in 
that  they  assume  residual  effects  even  during  the  first  period. 

Our  first  model  Is  as  follows:  an  observation  Y^^j,  taken  In  row 
1 and  column  J of  d,  can  be  expressed  as 

'‘l]  - “l+Bj  + ^(J(1,J)  + "ad-l.J)  + 'lj5  1 i 1 i P.  1 i J i n (1-1) 

where  are  uncorrelated  raindom  variables  with  common  variance 

and  zero  means.  The  unknown  constants  0^,  and 

Pd(i.i,j)  respectively  called  the  1^**  period  effect, 

experimental  unit  effect,  d(l,j)^*'  direct  treatment  effect  and 
d(l-l,j)^^  (first  order)  residual  treatment  effect.  We  assume 
that  for  1 ■ 1 ■ Pq(p,j)'  t**®  residuals  in  the 

first  period  are  Incurred  from  the  last  period.  This  way  each 
period  contains  residuals  from  the  Immediately  preceding  period. 


Mtualng  for  that  Matter  that  the  last  period  precedes  the  first. 

With  this  assumption  one  is  inclined  to  think  of  the  coluans 
of  d being  circular  and  hence  also  call  these  designs  circular. 
Another  Model  that  we  consider  is  the  same  as  (1.1)  but  no  period 
effect  is  present,  l.e., 

■*"  ^d(l,j)  **■  ^d(l-l,J)  ®ij*  (1*2) 

SiMllarly,  the  Model  May  not  contain  unit  effects, 

^ij  " ®i  ^d(i,j)  **■  Pd(i-l,j)  ■*“  (1«3) 

When  neither  period  nor  unit  effects  are  present  the  model 


assumed  is: 


^IJ  “ ^d(l,j)  ^d(i-l,j) 


(1.4) 


In  all  these  models  the  pn  observations  are  taken  on  a pxn 
array  d in  which  the  columns  are  assumed  circular,  as  described 
above.  Residuals,  being  the  principal  effects,  are  assumed  in  all 
the  models.  Each  of  the  last  three  models  can  very  often  be  associ- 
ated with  a sequence  of  observations  on  the  same  unit  (or  period) 
rather  than  a two  dimensional  array.  The  case  with  n>l  can  be 
particularly  meaningful.  When  n > 1 resting  periods  for  the 
unit  are  allowed  after  each  set  of  p treatments.  More  details 
on  the  Interpretation  of  the  above  models  can  be  found  in  Magda 
(1979). 

Let  us  say  a few  words  about  the  circular  nature  of  the 
columns.  One  can  certainly  encounter  practical  situations  where 


i : 
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I blocKs  are  of  auch  nature  that  residual  effects  are  transnltted 

[ along  one  direction  In  the  block  In  a circular  manner.  But  even 

If  the  natural  setting  of  the  experiment  Is  such  that  no  residuals 
would  appear  during  the  first  period  (as  Is  assumed  In  most  papers 

^ on  the  subject)  one  can  correct  that  In  the  following  way.  Intro- 

1 

I duce  a pre-period  (or  period  0)  In  which  treatments  are  given 

to  the  units  only  to  be  used  for  residuals  during  period  1.  Ho 
observations  need  to  be  taken  during  period  0.  But  It  is  essential 
that  during  the  last  period  of  the  experiment  the  saM  distribution 
of  treatments  as  In  period  0 will  be  made  to  the  units.  This  way 
the  residuals  from  the  last  period  (being  the  same  as  those  In 
period  0)  will  be  observed  during  period  1. 

Balancing  the  experiment  simultaneously  for  residuals,  periods 
and  units  becomes  very  natural  under  this  model. 

The  concept  of  strongly  balanced  uniform  RMD  has  recently 
been  brought  up  by  Cheng  and  Wu  (1979)  In  the  model  with  no  residual 
effects  during  the  first  period.  Necessary  divisibility  conditions 
for  the  existence  of  such  a design  are  (for  1^,  Xg  and  X posi- 
tive Integers ) : 


p ■ Xj^t,  n ■ Xgt  and  X • . 

O 

These  three  conditions  Imply  t |n  which  means  that  the  number 
of  experimental  units  can  get  rather  large.  In  the  case  of  a 
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p ■ and  n ■ 

turn  out  to  b«  the  only  neceaeary  divisibility  conditions. 


2.  fldf  initlOM » 

Let  us  recall  the  usual  notations  and  definitions. 

A repeated  nteasurenents  design  based  on  t treatments « n 
experimental  units  and  p periods  will  be  abbreviated  by  RND(t«n,p). 
The  collection  of  all  such  designs  is  denoted  by  04..^. 

tin»p 

Definition  2 .1.  A design  is  called  uniform  on  periods  if  each 
treatment  occurs  the  same  number  of  times,  say  in  each  period. 

A necessary  condition  for  this  to  happen  is  n « x^^t. 

Definition  2 .2 . A design  Is  said  to  be  uniform  on  units  if  each 
treatment  is  assigned  the  same  number  of  times,  say  X2«  to  each 
unit. 

This  can  occur  only  if  p > Xgt. 

Definition  2.^.  A design  is  called  uniform  if  it  is  uniform  on 
both  periods  and  units. 

Definition  ,2^.  A design  d € n«.  _ _ is  called  circular  strongly 
balanced  if  the  collection  of  ordered  pairs  (d(l,J),  d(l-fl,j)), 
l^i^p,  l^J^n  contains  each  ordered  pair  of  treatments 
(distinct  or  not)  the  same  number  of  tians,  say  X. 


i, 


i| 

1 


1 

1 


H 
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When  1 ■ p>  i-fl  ■ 1«  because  we  also  count  the  pairs 
(d(p»j)*  d(I,j)),  since  the  columns  are  assumed  circular. 

A circular  strongly  balanced  design  can  exist  only  if 
np  • \t^. 

It  is  not  hard  to  check  that  in  a circular  strongly  balanced 
uniform  RMD  the  following  divisibility  relations  hold: 

p « \^t,  n - ^ * ^1^2* 

Hence  X is  always  integral  if  X^  and  X2  are. 

Rxemnle  2 .1.  The  following  is  an  example  of  a circular  strongly 

balanced  uniform  RMD(4*4,8):  1234 

2 3 4 1 
4 12  3 
d*  - 3 4 12 

4 12  3 
2 3 4 1 

1234  (X  - 2). 

In  the  above  example  we  were  able  to  construct  a circular 
strongly  balanced  uniform  design  with  4 treatments « based  on  only 
4 units.  To  achieve  a similar  balance  on  the  4 units  when  no 
residual  effects  are  assumed  in  the  first  period  takes  at  least 
16  units.  (See  the  model  and  the  definition  of  a strongly 
balanced  uniform  RMD  in  Cheng  and  Wu  (1979)* 

Definition  2.^.  A design  d c 0*.  . „ is  called  circular  balanced 
if  the  collection  of  ordered  pairs  (d(i,j)«  d(i-i>l«j})»  1 ^ i ^ p* 


1 ^ J ^ n contains  each  ordered  pair  of  die tine t treatments  X 
times . 

Aa  before*  the  columns  of  d are  assumed  circular. 

The  divisibility  condition  for  the  existence  of  a circular 
balanced  RMD  Is:  pn  - Xt(t-l). 

1234567  1234567 

d-2345671  4567123 

4567123  2345671 

d e Is  a circular  balanced  RMD  uniform  on  periods. 


ijn  1 n ii'  1 1 M;i;iiir 


The  first  two  lemmas  which  we  shall  prove  have  easy*  but 
Important  statistical  Interpretations.  They  state  that  for  the 
purpose  of  estimating  linear  functions  of  certain  parameters* 
we  can  only  decrease  the  number  of  such  estimable  functions  and 
(In  connection  with  Lemma  2.1  of  Ehrenfeld  (1933))  only  decrease 
the  precision  of  our  estimators  by  allowing  more  parameters  In 
the  model. 

Under  the  general  linear  model 

E(y)  - +*“+Xj.9j.*  cov  Y - a^I 

we  denote  by  C^(9]^)  the  Information  matrix  of  G^* 
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For  the  two  linear  nodela 


E(Y) 


*■  ' n 

E and  E(Y)  . E I.e.,  cov  Y - a^I 

1«1  * 1«1  ^ ^ 


with  kin  we  have  0^(0^)  i Cj^Ce^). 


lEQfil:  For  convenience,  let  us  write  for  X^, 

1 i i i n,  and  reorder  the  parameters  In  the  design  matrix 
as  follows: 


. rxd)  x(i)  x<i) 

^1  * n »*n-l 


(1) 


If  0^  is  a vector  of  length  s,  the  sxs  matrix  In  the 
upper  left  hand  corner  of  Is  Cj^(0j) 

Do  the  following  operation  In  x(^)^X^^)  • (X(^)^x(l)j^ 

• l,n,n-l, , . . , 3,2 . Form  the  matrix  X^^)  x(2)  follows: 


I 


l( 
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x(2)'x(2)  . 


xU)'x(i)  x{i)'x(i)  xU)'x(i)  x<i>'x<i> 

11  1 n 1 n-1  •••  •*“3 


x(i)'xj^)  x<i)^c(i)  xCD'xd) 

n 1 n n n n-1 


X (i)'x(i) 
n 3 


xU)'xU)  xd)'x(i)  xd)'x(i) 

3 1 3 n 3 n-1 


x(i)'x(i) 
5 3 


Ixd)'x(i) 


x(i)Mi) 

ir«s‘> 


x(i)'x(i) 

3 2 


(X|1)';^1))-[^1)'41)  xi^)'xU)3 


where  (*1^^  Is  a generalized  inverse  of  x|^^  X^^^. 

Explicitly,  x(2)  . [4^^X^2)  x(2)^._^(2)^  x|2)].  where 


xp) , (i-4^h4^y 4^h-  4^y)4ii2^i^n^i 


and  ■ (I-(x|^)(4^^^xd))-  Then  In  the  same 


sxs  upper  left  hand  corner  we  have 
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C2(ei)  . x<2)'42)  . . 

. xU)'xU)(^i)'x<i)r  ^ xU)'x(i)  - C,(e,). 

We  can  now  proceed  recursively.  By  the  same  kind  of  elimina- 
tion In  X^**^  we  obtain: 


i • ^r^®l^  conclude  that  i Cj.(e^), 

for  all  1 i r i n-1.  Hence  i 1 i K i n. 


By  the  duality  which  exists  between  residual  and  direct  effects 
In  all  the  models  described  In  section  1,  it  is  clear  that  they 
both  share  the  same  Information  matrix  for  a fixed  model.  We 
denote  by  ^d^^  Information  matrix  for 

residual  effects  when  design  d is  used  under  models  (l.l)»  (1.2), 

(1.3)  and  (1.4).  Since  model  (l.i)  contains  more  parameters  then 

(1.4)  for  1 i 1 i 4,  Lemma  3.1  Implies 

i for  1 i 1 i 4 (3.1) 


We  shall  rely  on  this  fact  for  a unified  proof  of  several  results. 
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The  results  we  prove  concern  the  universal  optimality  of 
certain  circular  balanced  designs.  The  concept  of  universal 
optimality  has  been  Introduced  by  Kiefer  (1975).  It  Is  known 
that  If  a design  Is  \mlversally  optimal  then  it  Is  D-,  A-  and 
E-optlmal.  Using  the  sufficient  conditions  established  by  Kiefer 
(1975)  to  ensure  universal  optimality,  a design  d*  Is  universally 

i 

I optimal  for  residual  effects  In  0,.  _ _ under  model  (l.l)  If 

(1)  For  any  <1  € 0^  _ _ its  Information  matrix 
has  zero  row  sums. 


(ii) 


.(i) 


Is  completely  symmetric 


elements  of  C 


(1) 


are  equal  and  all  the 


are  equal. 


(lii) 


max  trci^^  i trC^P 

n n d 

t,n,p 


, l.e.,  all  the  diagonal 
off-diagonal  elements  of 


We  proceed  In  showing  that  (i),  (ii)  and  (ill)  are  satisfied 
for  suitable  choices  of  d*  under  the  various  models . 

In  order  to  show  that  has  row  sums  zero,  for  1 H 1 4, 

we  first  prove  a lemma.  This  lemma  can  be  also  derived  from  Lemma 


2.1  of  Ehrenfeld  (1955). 
lesm  Let  A,B  0. 

If  A ^ B then  the  row  span  of  A Is  contained  In  the  row 


span  of  B. 


-n- 


ProQf  1 We  can  write  A * N'H  and  B - M'M  where  N and  M 
are  of  full  row  rank.  Than 

MX  - 0 • x'M'Mx  - 0 « x'N'Nx  - 0 « Nx  - 0. 

Therefore  (x:Nx  ■ Oj  e (x:Nx  « Oj,  or  the  orthogonal  complement 
of  the  row  span  of  M la  contained  in  the  orthogonal  complement 
of  the  row  apan  of  N.  Hence 

(row  apan  of  M)  s (row  apan  of  N). 

Then  the  column  apan  N*  M ■ B containa  the  column  apan  of 
N'N  ■ A.  Since  A and  B are  aymmetric  the  aame  ia  true  for 
the  row  apana. 

Note  that  if  the  row  auma  of  are  zero,  then  the  row 

(4) 

apan  of  ' ia  orthogonal  to  the  vector  1,  whoae  entriea 

are  all  1.  Since  by  (5.1)  i for  1 i i i 4,  Lemma 

5.2  aaaurea  that  the  row  apan  of  ia  alao  orthogonal  to 

1 aince  it  ia  contained  in  the  row  apan  of  In  order  to 

ahow  (i),  for  1 ^ i ^ 4,  it  therefore  aufficea  to  ahow  that 

~ Ik) 

the  row  auma  of  ' are  zero. 

It  ia  eaay  to  eatabliah  that  “ ^d  “ **d®d^*^* 

®d  * ^**di^  **d  " ^®dij^  where 

^di  * number  of  appearancea  of  treatment  i in  d. 

®dij  * number  of  timea  treatment  i ia  immediately 

preceded  (columnwiae)  by  treatment  J in  d. 
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The  following  dependencies  hold  for  any  d e ^ 

t,n,p 


t 

T 

J-1 


^’dl  “ 4^,  “dlj  * “djl 


and  r r 
1-1 


dl  • P" 


The  (1,J)^^  entry  of  “ ®d  “ **d^d^*?i 


(3.2) 


" kfl  »’dK 


where  6^.  - 1 If  1 - J and  0 otherwise.  Indeed*  the  sum 

(4) 
d 


along  the  1^^  row  In  Is  (using  (5.2)): 


E (6.  .r.  - E - E E ®^^lr®H^b. 

-1  ^ K-1  *^dK  J-1  K-1  *^dk 


k^l  'dk 


“dlk  “djk 


'■<‘1^  ■ kfl  '^dK 


’ ’^dl  “ “dlk  * ‘‘dl  " ‘‘dl  • 

« 

In  what  follows*  d^^  denotes  a circular  strongly  balanced 
uniform  RMD*  ^2  denotes  a circular  strongly  balanced  uniform 
on  units  RMD*  d^  denotes  a circular  strongly  balanced  uniform 
on  periods  RMD  and  finally*  d^  denotes  Just  a circular 
strongly  balanced  RMD. 


rrs 
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By  doing  the  standard  row  operations  under  the  four  models 
we  establish  that 


(I  - J)  for  1 i 1 i 


(3.3) 


where  1 Is  the  Identity  matrix  of  order  t and  J is  the 
txt  matrix  with  all  Its  entries  equal  to  1.  Note,  for 
example,  that 


C<i)  . f (I  - i J) 


when  J ^ 1,  but  not  when  J > 1. 

The  condition  (11)  on  complete  symmetry  Is  contained  In 
(3.3).  It  remains  to  check  (111),  l.e.. 


,(i) 


deO 


max  trCi  ' ^ trC'i'  for  1 i 1 ^ 4. 


(i) 


t,n,p 


This  also  happens  to  be  easy  to  show,  by  making  use  of  both 
(3.1)  and  (3.3).  Prom  (3.1)  and  (3.2) 

2 


trc<l>  i trof")  - tr(D^  - - pn  - ^ , 

for  all  1 i 1 i 4.  Prom  (3.3) 

trC^*^  - pn  - for  all  1 i 1 i 4. 

A . V 


t m! 


It  therefore  suffices  to  show  that  d^  minimises  £ I . 

J-1  k-1  *^1 


dk 


a 
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over  Ot.n.p  1 i 1 i 4.  Indeed,  since 

t t 

jfl  “dJK  - 'dk  •"*  *‘dlc  • P" 


r t d.ik  i i 2 . 

E E -T^  - E E "Vnt  i 

J-1  k-1  *^dk  K-1  *^dk  J-1 


t 

E 

k-1 


k-1  ® 


sn 

t • 


The  inequality  used  above  corresponds  to  the  fact  that  the  sum 
of  squares  of  t nonzero  terms  having  constant  sum  Is  minimized 
when  they  are  all  equal. 

Having  satisfied  the  conditions  (1),  (11),  and  (ill)  for 
1 ik  1 4 we  have  proved  the  following  theorem: 


Theorem  ^.1.  Under  the  model  (1.1)  whenever  a design  d][  exists. 
It  Is  universally  optimal  for  the  estimation  of  direct  as  well 
as  residual  effects  over  the  collection  of  designs  with  the  same 
parameters  (1^  4). 

The  following  stronger  optimality  result  holds  If  we  suitably 
restrict  the  class  of  competing  designs: 


ThgQraa  ^i2«  under  the  model  (1.1)  whenever  a design  d^  exists. 
It  minimizes  the  variance  of  the  best  linear  unbiased  estimator 
of  any  contrast  of  direct  effects  and  any  contrast  of  residual 


I 


1 


{ 
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effectt  over  the  collection  of  equlrepllcated  deeigne  with 
the  seme  psrsmeters  (1  ^ 1 4). 

Proof!  A design  d c called  equlrepllcated  If 

I,  or  equivalently,  if  m ,,,  m Let 

« 

^t,n,p  the  collection  of  equlrepllcated  designs  with 

parameters  t,n  and  p.  By  Lemma  2.1  of  Ehrenfeld  (1953),  It 
la  enough  to  show  that 


^ A* 

for  all  d c 0,.  _ _ and  all  txl  vectors  x.  We  can  write 
X - y + 1,  where  y'l  - 0.  Then  knowing  that  Is  sym- 

metric and  cl^^l  - 0 for  all  d € oil  ^ we  have 
d mt  t,n,p 

x'c^^^x  - (y'  + (y  + 1)  - y*C^^^y 

where  jr'  1 • 0.  It  suffices  therefore  to  show  that 

m 

y#C^i)y  ^ y^C^^^y 

^1 

for  all  y such  that  y^l  ■ 0.  We  know  from  (3.1)  that 

i.  cj"*  . D^,  - i Dj  . I 


for  d c Hence 

/C^^^y  ^ r ^ 1 y • y*  ^(I-'Jj)y-  y'c^jV. 

d. 
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•ince  y'  Jy  ■ 0.  This  proves  the  theorem. 

Up  to  this  point  we  have  seen  that  under  various  models 

various  kinds  of  circular  strongly  balanced  designs  proved  to 

be  optimal.  We  now  turn  to  the  concept  of  a circular  balanced 

RMD  and  show  soise  of  its  optimal  properties. 

« 

Let  denote  a circular  balanced  uniform  RMD  (see 
Definition  (2.5)*  ^3  • circular  balanced  uniform  on  units  RMD, 
5*  a circular  balanced  uniform  on  periods  RMD  and  finally, 
ej  Just  a circular  balanced  RMD.  As  before,  will  denote 

the  information  matrix  for  residual  (and  hence  also  direct) 
treatment  effects  when  model  (l.i)  is  assumed.  As  previously 
we  have 

i for  liii4  (3.4) 

with  ■ D^  - where  D^  and  M^  are  as  before 

and  the  relations  (5.2)  hold.  One  shows  that  has  row 

U 

suns  zero  either  directly  or  using  Lemma  5.2  as  we  did  previously. 
By  carrying  out  the  standard  row  operations  it  turns  out  that 

■ ("  • - 1 ’’)  1 i i 

It  is  now  convenient  to  use  the  following  Lemma  in  order  to 
establish  trace  maxinality.  We  Just  state  it  and  omit  the  easy 
proof. 
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2 


i.  Let  be  n Integers.  The  minimum  of  £ mf 
n 

subject  to  I m.  • k is  sttsined  when  k - n intf-^^  of  the 
1-1  ^ 

m^  equal  1 + n - k + n lnt(^  of  the  m^^  equal 

int(j[)*  corresponding  minimal  value  is 


•"(‘"‘(is))  + <*'‘-">(‘"‘(S)  + 


The  symbol  int  x symbolises  the  integral  part  of  the  real 
number  x. 


Turning  to  the  maximisation  of  the  trace*  by  (5.4)  and  (5.2) 


t t m: 


it  suffices  to  show  that  6j.  minimises  E E . In  this 


■'ll  ww 

^ J-1  k-1  *^dk 

case  a solution  with  all  equal  and  integral  is  not 

possible  because  of  the  present  divisibility  condition.  Hence 

we  have  to  carry  out  the  minimisation  over  Integers*  using  Lemma 

3.5*  l.e.* 

'^JK  i -‘(‘"‘(^))^  + - t)(lnt(^))  + r, 

.r. 


"dk 


i -(‘"‘("t‘))<'’<lk  - *>  + "di.  ^ "dk- 


For  any  <1  c 0*.  _ . with  X - 1 we  therefore  have: 
b *n*p 


t m' 


41r  4 


£ ■: 


_ — «•  _ tm  av*  4W  ^ £ W 

J-1  k-1  *^dk  k-1  Mk  J-1  k-1 


dk 


"dk  i 


with  equality  for  (with  X - 1).  This  establishes  the 

universal  optimality  of  bt  over  0*.  _ - when  pn  - t(t-l) 

•*  b*n*p 

(or  equivalently*  when  X - 1). 


t 

n 


1 i 


J 


(■ 

AM 


In  eaa«  X > 1*  C c it  ndnimlted  by  &•.  only  in  a 

J-1  k-1  *^dk  ^ 

aaaller  clatt  of  dttigna.  Let  Ql  _ » the  collection  of 

vf  HfP 

detignt  d c ()t«n«p  ^'^th  the  property  that  no  treatment  pre- 
cedes itself,  i.e.,  m ^ ^ ■ (d:  - 0,  1 i i i tj. 


Then  for  any  d c nl 


with  equality  for  t 


The  optimality  of  0^»l^i^X,  cannot  be  extended  to  all 
n n X > 1*  The  following  example  shows  this. 


6^  is  a circular  balanced  uniform  on  periods  RND  with  X ■ 5, 
while  dt  is  a circular  strongly  balanced  uniform  on  periods  RND 
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The  both  belong  to  By  Theorem  3.2,  d;  Is  better  than 

63.  This  shOMs  that  the  optimality  of  6*  can  only  be  proved 
over  a attbclftag  of  0^  p when  1 > 1. 

Let  us  now  state  the  results  established  above. 

TheorAffl  ^.3.  Under  the  model  (1.1)  whenever  a design  6^  with 
X - 1 exists,  it  is  universally  optimal  for  the  estimation  of 
direct  as  well  as  residual  effects  over  the  collection  of  de- 
signs with  the  same  parameters  (1  i 1 ^ 4), 

XhsJlIclLiJ.  Wh«never  a design  6*  e exists.  It  Is  uni- 

versally optlaal  for  the  estloatlon  of  direct  as  well  as  residual 
effects  under  model  (l.i)  over  0^j^p(lil^4), 
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